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Abstract 

Within the framework of generalized combinatorial approach, the complexity is 
determined for infinite set of self-similar hierarchical ensembles. This complexity 
is shown to increase with strengthening of the hierarchy coupling to the value, 
which decreases with growth of both scattering of this coupling and non-extensivity 
parameter. 
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Hierarchical structure is one of the most universal peculiarities of complex 
systems in physics, biology, economics and so on [l]-[5]. Their evolution is 
shown [6] to reduce to anomalous diffusion process in ultrametric space of the 
self-similar hierarchical system, whose steady-state distribution over hierar- 
chical levels is given by the Tsallis power-law, being inherent in non-extensive 
systems [7] . Well-known feature of hierarchical systems consists in that every 
of statistical ensembles of given level breaks with passage into lower one to 
subensembles, then every of these breaks to more small subensembles of the 
following level, and so on (see, for example, Figure 1). From the statistical 
point of view, the set of above (sub)ensembles is characterized by the com- 
plexity, whose value determines the scattering of the hierarchical coupling - 
in analogy with the entropy in usual statistical systems. This article aims to 
determine the complexity for self-similar hierarchical ensembles. 

According to Ref. [7] the non-extensive statistics is based on the definitions of 
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Fig. 1. Characteristic form of the internal energy relief (a) and related hierarchical 
tree (b) of complex system [3] 

both logarithmic and exponential functions by the equations 



x — 1 

\n q (x) : = — - , exp g (a;) := [1 + (1 - q)x]+ q ; [y} + : = max(0,y), q < 1, 



that are reduced to the usual functions in the limit q — > 1. Introducing g- 
deformed production and ratio of positive values x, y as follows: 



x® q y 



i 

1-9 



x ^-y l -i + !}]-"; x,y>0, (2) 



it is easily to convince that these satisfy to usual properties lia q (x £g> g x) = 
ln q x + ln q y, ln q (x Q) q x) = ln q x - ln q y; exp g (a;) <S> q exp q (y) = exp q (x + y), 
exp q (x) <Z) q exp q (y) = exp q (x - y). 

Within the framework of combinatorial approach [8] , the g-deformed statistics 
is reduced to consideration of the generalized factorial N\ q := I ® q ■ ■ ■ ® q N , 
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whose logarithm is as follows: 

\n q (N\ 
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In thermodynamic limit A" — > oo, above sum is estimated by related integral, 
whose calculation gives 



£-\n q N-£- q + 0(\n q N), q^2, 
N-\nN+0(l), q = 2. 



(4) 



Defining g-deformed multinomial coefficient by the equation 

= (N\ q ) (d q \W. q ) (JV*!,)], 
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where the set of integers Aj satisfies to the condition A" = Y%=i Ni, we find 



' N 



N JVi 
1=1 ii=l 



iVfe 



E 



1-9 



ifc=l 



1/(1-9) 



From here, similarly to Eq.(4), one obtains the expression 
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for the Tsallis entropy 

Cg(pi,...,pjv) := -EPi ln 9Pi 
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Above formalism is generalized easily for consideration of the hierarchical sys- 
tems [9]. Assume that the A" states of upper level are distributed over the 
ensembles % = 1, . . . ,n, every of which contains Aj states. In turn, Aj states 
are bunched in m 8 subensembles ij, every of which contains Nij states, where 
the relations £^i N tj = N h E?=i N t = N are fulfilled. Then, instead of the 
multinomial coefficient (5) it is necessary to use the expression 



N 

N u ...N n 



< N ^ 



iVi 



N„ 



(9) 



3 



whose g-logarithm is 

s U . "O, = s U "J. + I s U .\J,- (io) 

As a result, using estimation (7) arrives at the connection between the com- 
plexities of the nearest hierarchical levels 

C Q (pn, • • • ,Pnm n ) = Cq (pi, . . . ,p n ) + ^2pfC Q (—,..., — I , (11) 

i=i VP* Pi / 

whose distributions over states are given by the equations pij = Nij/N, pi = 
Ni/N (here we introduce 'physical' non-extensivity parameter Q = 2 — q, 
1 < Q < 2). Using definition (8) and connections pi = Y,™LiPiji at condition 
Pi — Pij t <C Pi gives estimation 

C Q {^...M^t{ P -^)\ (12) 

V Pi Pi ) 2 j% V Pi J 

with whose accounting one finds 



Cq (pn, • • • ,P„ m J - C Q (pi, . . . ,p n ) « — ^p^ S (Pi - Pi*) • ( 13 ) 

=i ii=i 
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If statistical states are distributed within microcanonical ensembles, then both 
probabilities and related complexities are determined by the level number n: 

{Pi]l => Pn, {Pin)T => Pn+U Cq (pi, ...,p n ) => C(n), Cq (pn, . . . ,p nm J =>• 

C(n + 1). As a result, Eq.(13) takes the simplest form: 

C(n + 1) - C7(n) « ^y^p?" 2 (Pn - Pn + i) 2 . (14) 



At fixed level number n 1, relation obtained is reduced to the differential 
equation 

<9 2 C(n) 



whose integration gives the dependence 



= Q(n+l)p-^\ (15) 



C(n)=rf-^_p?- 1 Pfl + ^±ipJ (16) 



at the boundary conditions 

9C(ra) 



C(n = 0) = 0, 



<9p r , 



= 0. (17) 

n=0 
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Then, using self-similar distribution [6] 

Pn = A [A + (Q - l)(n + 1)]-^ , A = (2 - Q) [(Q - 1) + A]3f* , (i 8 ) 

where parameter A determines the scattering over hierarchical levels, we ob- 
tain the complexity distribution shown in Figure 2. As it should be, on the 
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Fig. 2. Complexity dependence on hierarchy level number (curves 1, 2, 3 relate to 
parameters A = 0.1, 1.0, 10; Q = 1.5) 

upper level n — 0, where only one statistical ensemble takes place, the com- 
plexity is C = 0. With passage to lower levels the value C increases, taking 
maximum value at no = (Q — 1) + A, after which characteristic magnitude 
Coo = Po is reached. However, it follows to take into account that above con- 
sideration is applicable for level numbers n ^> 1 only. Therefore, the pointed 
out maximum is displayed at condition A ^> 1, when distribution (18) differs 
very weakly from the exponential one, being inherent in usual additive sys- 
tems. At moderate scattering of the hierarchical coupling A ~ 1, the system 
complexity increases fast from zero to the limit value 



C n 



2-Q 



(Q-l) + A 



(19) 



which decreases with growth of both scattering A and non-extensivity param- 
eter 1 < Q < 2. 

Evolution process of the hierarchical ensemble is determined by both diffusion 
time Td = (A 2_< ^ / 'D)rfi W and scale r = n 2 T of the steady-state reaching, 
where D is diffusion coefficient, r is time of passage between nearest levels 
[6]. At t r d , the diffusion over the hierarchical tree has anomalous character 
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to be determined by the law n(t) ~ Q 1 ^(t/r ) 1 ' ,( ^. With time growth within 
the interval ~ t <C r, when contributions of anomalous drift and diffusion 



are comparable, the passage into the normal regime n{i) ~ y 2(t/Vo) happens. 
At t <C r, the distribution over hierarchical levels tends to the steady-state 
law (f8) according to the time dependence 

Vait) = Pn [l-(t/T n )-Q^] , (20) 

which is reduced to the value p n (t) = at t < r n , r n = (A/QD)nro, n ^ 
(here, initial distribution p n (t = 0) = 5 n0 is taken in form of the Kronecker 
5-symbol). 

Non-stationary complexity is determined by Eq.(15), where one should take 
the boundary conditions 

C(n, t = 0) = 0, C(n, t = oo) = C(n), n ^ (21) 

instead of Eqs.(17). As a result, we derive the time-dependent complexity in 
the form 

C ^ t )=(f i -^ 1 Po' 1 )pn(t) + ^\pm, n»l, (22) 

being generalization of the steady-state distribution (16). Inserting Eqs.(18), 
(20) into Eq.(22), we obtain the time dependence of the complexity shown in 
Figure 3. During the ballistic interval t < r n , the complexity keeps the initial 
value C — 0, whereas in the course of the time t > r n it increases fast before 
the steady-state value (16), being bounded by the maximum value (19). 
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Fig. 3. Time dependencies of the complexity on different hierarchy levels (their 
numbers are noticed near curves) at parameters A = 1.0, Q = 1.5, D = 1 
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